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, Abstract. We review some recent progress in quantum field theory in non-commutative 

space, focusing onto the fuzzy sphere as a non-perturbative regularisation scheme. We first 
introduce the basic formalism, and discuss the limits corresponding to different commutative 
or non-commutative spaces. We present some of the theories which have been investigated 
, in this framework, with a particular attention to the scalar model. Then we comment on the 

1^ ' results recently obtained from Monte Carlo simulations, and show a preview of new numerical 

data, which are consistent with the expected transition between two phases characterised 
by the topology of the support of a matrix eigenvalue distribution. 
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1 Introduction 

^ ■ 

The study of quantum field theory in non-commutative spaces [H EJ [Sj H] is based on strong 
physical motivations, since this formalism is relevant to the description of the quantum Hall 
effect [5], to certain aspects of string theory [H O [HI [9l [TOl [11] , as well as to a possible formulation 
5^ I for a quantum theory of gravity |12| . 

The most extensively studied non-commutative space is probably the Groenewold-Moyal 
Mg space, defined as the algebra of functions Ag{W^) generated by the coordinates (/i e 
{1, 2, . . . , n}) with the commutation relations: 

[x^x'^], = ^r^ (1) 

where, for simplicity, 6^'^ can be taken to be independent of x, and * denotes the Groenewold- 
Moyal star-product. For functions /, g in Ae{W^), the latter is defined to act as [13 ^ 114 1 [TS]: 

f{x)*g{x) = f{x)e^2^-'"'^^g{x). 

A number of dramatic implications stem from equation ([1]), which have been studied in detail in 
the literature [I6l[171[l8l[ia[20l[2ll[l5l[22l|23l[a[25l[2^ In particular, 

some issues which are of strong relevance from the physical point of view include: 



*This paper is a contribution to the Proceedings of the O'Raifeartaigh Symposium on Non-Perturbative and 
Symmetry Methods in Field Theory (June 22-24, 2006, Budapest, Hungary). The full collection is available at 
|http://www.emis.de/journals/SIGMA/LQR2006.html| 
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• The "ultra- violet /infra-red (UV/IR) mixing" phenomenon: The effective action describing 
QFT in a Groenewold-Moyal space diverges when the components of the external momen- 
tum along the non-commutative directions vanish. This effect arises from the integration 
of high-energy modes in non-planar loop diagrams, and is related to the non-locality prop- 
erties of quantum field theory in a non-commutative space. 

• Renormalisability: One of the potentially most dangerous consequences of the UV/IR 
mixing phenomenon is the fact that it threatens the possibility to renormalise quantum 
field theories defined in a non-commutative setting. However, the proof of perturbative 
renormalisability at all orders has been eventually worked out adding an oscillator term 
to the potential of the scalar model [281 [29l [30| [3T1 [32] . 

The theoretical predictions for these models can be compared with non-perturbative results 
obtained from numerical simulations. In fact, it is possible to define a well-suited lattice formula- 
tion of the non-commutative theory as a unitary matrix model [33], which possesses a continuum 
limit characterised by finite space volume and finite non-commutativity. This formulation al- 
lowed to investigate numerically various aspects of these models [3 ^ [35l[36l[37t [38 ] [39l BOl HT] H2] . 

A different class of non-commutative spaces is given by fuzzy spaces; the basic idea underlying 
their construction is to use a finite-dimensional matrix algebra to approximate the infinite- 
dimensional algebra of functions on a manifold. This can be done on even-dimensional co-adjoint 
orbits of Lie groups, which are symplectic manifolds [I3l[lllll5lll6l[171ll8l[l9l[50l[51], like the 
two-sphere 5^ and the CP" complex projective spaces. 

In the following discussion, we shall concentrate our attention onto the fuzzy two-sphere 
Sp [13]; it provides the simplest example of a fuzzy space, and depends on two parameters: the 
dimensionful, real and positive, radius R, and a dimensionless, integer- valued cut-off, which is 
associated with the value of the maximum angular momentum ^max- In particular, the ordinary, 
commutative algebra of functions on the sphere is replaced by a non-commutative counterpart, 
obtained truncating the expansion of functions in the basis of irreducible representations of the 
su{2) algebra to a maximum angular momentum /max- This geometrical construction of Sp 
offers a natural representation for functions on the sphere in terms of matrices belonging to 
MatAr(C), with N = Z^ax + fl 

The fuzzy approach explicitly maintains the symmetries of the original manifold for every 
value of A^, and it accommodates various topological features in a natural way - see [50] and 
references therein for details. It provides a well-defined, non-perturbative regularisation scheme 
for quantum field theory, and - at least in principle - can be considered as a potential alternative 
to the lattice discretisation. The representation of fields on the fuzzy sphere in terms of finite- 
dimensional matrices also opens up the possibility to study interacting theories via numerical 
simulations: this has motivated various recent works [521 ESI EH ESI IMl ET] [58] [59]. 

Indeed, the fuzzy models reveal a rich and complex structure, encoding a number of interesting 
properties, which show up even for the simplest theories - e.g.: for a scalar field with quartic 
interactions. 

The aim of this paper consists in presenting an overview on some of these aspects, comparing 
theoretical expectations and results of Monte Carlo simulations, and highlighting their physical 
interpretation in different limits. 

The sections below are organised according to the following structure. In Section [21 we 
first introduce the general formalism for the construction of the fuzzy sphere, and discuss its 
different geometrical limits; then we define the description of QFT on the fuzzy sphere, and 
present the theoretical predictions which are expected in different limits. Next, the theoretical 
predictions are compared with non-perturbative results obtained from numerical simulations in 

^This construction can be easily obtained by taking the tensor product of two vector representations of angular 
momentum ^^j^. 
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Section [3l Finally, in Section U] we discuss the relevance of these results for the program to use 
the the fuzzy approach as a regularisation scheme for quantum field theory in commutative or 
in non-commutative spaces, and conclude with some brief remarks. 



2 General formalism and theoretical predictions 

The construction of a fuzzy space is based on the quantisation of the algebra of functions 
defined on the corresponding (commutative) manifold. For the two-sphere, the construction 
can be summarised as it follows [60]: The infinite-dimensional algebra of polynomials in the 
{xi}i=i^2,3 commutative coordinates on the two-dimensional sphere XiXi = E? embedded in 
is replaced by the non-commutative algebra generated by {xj}j=i,2,3 operators, which obey the 
(rescaled) commutation relations of the su{2) algebra: 

[xi,Xj\ = -===ieijkXk, with: ^^x1=R^. (2) 

i=l 

On the fuzzy sphere, the counterpart of functions defined on the commutative sphere are 
matrices belonging to MatAr(C); a natural mapping can be defined truncating the basis of 
harmonic functions on 5^ to a maximum value for the angular momentum Zmax, and associating 
the harmonic functions with total and third-component angular momentum quantum numbers I 
and m to the {^i,m} polarisation tensors satisfying: 



[Li, =l{l + [L^Am] = m% 

where the {ii}i=i,2,3 are the generators of the su{2) algebra in the N x N matrix representation. 
This mapping choice, however, is not unique - a convenient alternative mapping being the one 
defined in terms of coherent states. 

Integrals of functions on the commutative sphere are replaced by the matrix trace operation: 
given any <I> and ^' elements in MatAr(C), an inner product on Sp is built as: 

($,M/) = ^tr(<I>tM/). 

The vectors describing derivations on the commutative sphere are represented through the 
adjoint action of the {ij}j=i,2,3 generators: 

It is important to discuss what is the meaning of the "commutative limit" in this context: 
The matrix algebra underlying the fuzzy space is truly non-commutative for any N; however, 
the right-hand side of equation ([2]) vanishes in the — > oo limit, keeping R fixed. Therefore, 
given a smooth function defined on the commutative sphere, whose expansion in the spherical 
harmonics' basis involves non-negligible coefficients up to a certain angular momentum only, 
the effect of this discretisation will be vanishing in the large- limit at fixed R. 

So far, the interpretation of this limit is only concerned with the geometric properties of the 
manifold - not with the physical properties of a quantum field theory defined in this setting. 

Indeed, the quantum features of a physical field theory spoil the correspondence between the 
fixed-i2, large-A^ limit of the (naive formulation of the) theory on the fuzzy sphere, and the 
corresponding model on the commutative manifold [61] . 

This can be seen in a perturbative calculation for the (real) scalar model with quartic inter- 
actions; for this model, the Euclidean action on the fuzzy sphere can be written in the following 
form: 

S ($) = f tr (i [L„ <&]t [L„ + ^ci>2 + |$4^ (3) 
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(where $ G Mat7v(C) is taken to be Hermitian). At one-loop order in a perturbative expansion, 
the non-commutative nature of the model manifests itself in the comparison between the contri- 
butions to the one-particle- irreducible two-point function from planar (/p) and non-planar (/np) 
diagrams - they differ by a quantity which is finite, non-vanishing, and which has a smooth 
dependence on 



N- 



/np — -^p 



ij{j + l){2j + l 



When the N oo, fixed-ii, limit is taken, this difference between I^p and Ip has an impact on 
the one-loop effective action, as it induces a weakly non-local, /-dependent, deformation of the 
dispersion relation on the fuzzy sphere. This effect has no analogue in the commutative setting, 
and is therefore called the "non-commutative anomaly". 

In order to recover the correct model in the commutative limit, one should redefine the inter- 
action term in in equation ([3]) via a normal-ordering prescription [62], cancelling the undesired 
momentum-dependent quadratic terms in the effective action; alternatively, one could include 
rotationally symmetric higher derivative terms in the matrix action [63]. 

On the contrary, studying the model defined by the (naive) matrix action in equation ^ , one 
expects to observe deviations with respect to the quantum theory on the commutative sphere. 
Although the derivation of the non-commutative anomaly in [61] was worked out for the /x^ > 
case, one may expect that a similar phenomenon occurs when ij? < 0, too. 

For QFT on the commutative sphere, the latter regime is particularly interesting, as it corre- 
sponds to a situation in which, classically, the $ ^ — # symmetry of the model is spontaneously 
broken. The situation from the quantum point of view is different, however, because - due to 
the finite size of the system - tunneling events characterised by a finite Euclidean action connect 
the classical vacua, and the true quantum ground state is actually unique for any finite value 
of the parameters. Therefore, strictly speaking, one does not expect phase transitions on the 
commutative sphere with fixed, finite radius; this, however, might no longer be the case for the 

— > oo limit of the fuzzy model, if "anomalous" effects take place. 

For the quantum theory on the commutative sphere, it is also interesting to note that the tun- 
neling phenomenon between the (uniform) classical vacua can be mediated by typical quantum 
events which break the SO (3) symmetry - for instance, through a p-wave field configuration 
(with axial symmetry only). A signature of these events would show up through non- vanishing 
expectation values for the modes associated with non-vanishing momenta. In the commutative 
setting, this phenomenon is possible because of the finite size of the system; on the contrary, one 
would not expect a non-uniform breaking symmetry pattern in the infinite commutative plane. 

Another important implication of equation is that the contribution of non-perturbative 
diagrams to the two-point function is not divergent for vanishing values of the momentum. This 
has to be compared and contrasted with the non-commutative plane setting, where the UV/IR 
mixing phenomenon occurs. Yet, the divergence which is encountered in the the Groenewold- 
Moyal space is exactly recovered, once a different scaling limit of the fuzzy sphere is considered: 
this is discussed below. 

A local approximation to the non-commutative plane can be obtained, in the following 
double scaling limit of the fuzzy sphere: 

oo, R oo, with: 9 = — ^ fixed. (5) 

This can be seen considering a stereographic projection map from the sphere to the (2/1,1/2) 
plane: 

y+=yi+ iy2 = 2Rx+ {R - £3)'^ , y- = yi - iy2 = 2R{R- £3)'^ x_. 
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In the hmit iji and y2 satisfy: 

and the non-commutative anomaly induces a logarithmic infrared divergence proportional to 
log {p9) in the non-planar contribution to the IPX two-point function. This divergence is exactly 
the same which is found in the perturbative calculation for the quantum theory defined in the 
Groenewold-Moyal plane [22], therefore in this case the (double scaling) large- A'^ limit of the 
fuzzy sphere reproduces the expected result, and the fuzzy space can be used as a regularisa- 
tion for non-commutative spaces. In principle, the construction can be generalised to higher 
dimensional cases, using fuzzy CP" spaces, or direct products of fuzzy spheres. 

The regularisation of Euclidean scalar field theory in even-dimensional, non-commutative 
spaces by means of fuzzy manifolds was used in [63] ■ There, the properties of the model 
were investigated non-perturbatively in a matrix representation for the field, and a critical line 
corresponding to a phase transition in the associated eigenvalue distribution was derived in the 
large-matrix limit. 

This transition is a change in the topology of the support of the eigenvalue distribution: for 
values of //^ above a critical value ("single-cut phase"), the eigenvalue distribution p{(j)) has 
a connected support of finite width, centred in zero, and would reduce to Wigner's semi-circle 
law in the non-interacting case: 



1 

2^ 



1 



5off ( 0" + - ) + /^L 



yr^, with: = ^ (4 - /"D , (6) 



where g^s and /U^^ depend on (the renormalised values of) g and /i^. 

By contrast, for values of /i^ which are negative and large in modulus, there exist two sym- 
metric, disconnected regions of the real axis with non- vanishing probability for the eigenvalues 
("two-cut phase"). In particular, the probability density for the zero eigenvalue vanishes when 
/i^ is decreased down to the critical point, and for lower /i^ values a gap opens up in the dis- 
tribution. This phenomenon (which would not take place for a model in ordinary commutative 
space) is related to the UV/IR mixing effect; it is predicted on the basis of a matrix model 
description, which relies on the dominance of planar diagrams over the non-planar ones, and 
can be interpreted as a manifestation of the "striped phase" [24^ I25j. 

Figs. [1] to m show p{(j)) for different cases. Fig. [1] corresponds to the pure Wigner's semi-circle 
law; next. Fig. [2] displays the eigenvalue density profile in the "one-cut" regime. When /x^^ is 
decreased down to a critical value, p{Q) vanishes - Fig. [3]- and for lower values of the support 
of the distribution breaks down into two disconnected components (Fig. [4|). 

The arguments underlying the derivation in [6l] are expected to be stronger in the four- 
dimensional case than in the two-dimensional setting, because iw D = 2 the dominance of the 
planar diagrams over the non-planar ones is weaker; however, even in D = 2, one can still look 
for a signature of the expected effects - and, as it will be discussed in Section [3l the agreement 
with the numerical results appears quite remarkable. 

A different model was studied in [5l]: it is a scalar field theory in three-dimensional space, 
in which two coordinates are regularised using a fuzzy sphere, while the third one (labelled as t) 
is considered as purely commutative, and is regularised via a conventional lattice discretisation: 



1 

+ 2 



^ ^ + — <I>2(t) + -^\t) \ . (7) 
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Figure 1. Eigenvalue distribution correspond- Figure 2. Typical profile of the eigenvalue distri- 

ing to Wigner's semi-circle law: this is obtained bution in the "one-cut phase" . 

setting /i^ff = 4 in equation ([6]). 

Also in this case, different limits can be investigated, and one may expect a general pattern 
qualitatively similar to the lower-dimensional case discussed above. 

On the other hand, a fuzzy space description for gauge theory was studied in some papers, 

including [BSllMlEZlEHllMlllQlEIlliaEsllHIZa^ 

Supersymmetric versions of these models have been proposed, too. In particular, a supersym- 
metric generalisation of the fuzzy sphere was discussed in [771 EH CSl ED] , essentially based on 
the osp{2\l) algebra. Other works where the interplay between supersymmetry and the features 
of fuzzy models were discussed include [8ll[82l[83l[5Sl[7Il[Ha[85l[86]. 



3 Focus on the scalar model in two dimensions: 
a comparison with numerical results 

The models described in Section [2] can be quantised in the path integral approach, defining the 
expectation values of generic observables O = 0{ci^m) as: 

where the ci^m are the components of the ^ matrix in a given basis (e.g.: in the polarisation 
tensor one), and can be interpreted as the dynamical degrees of freedom of the matrix model. 

The right-hand side of equation ([8]) can be evaluated perturbatively, or estimated numerically 
from Monte Carlo simulations, in an approach similar to standard lattice field theory computa- 
tions, namely: via averages over a sample of (independent) configurations with statistical weight 
proportional to exp [—S ($)]. 

In this section, we summarise the present status of numerical results for the models discussed 
above, with a particular attention to the scalar theory on the two-dimensional fuzzy sphere. Our 
aim is to discuss the general physical interpretation of these results; therefore, we encourage the 
readers to refer to the papers mentioned, for precise details about the results, and for most of 
the technical aspects associated with the simulations and with the measurement of different 
observables. 

The real scalar model with quartic interactions on the fuzzy sphere was studied numerically 
first in [52], where it was pointed out that the model exhibits three different phases: a disorder 
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Figure 3. At the transition point, the probabih- 
ty to observe the zero eigenvahie becomes vani- 
shing. 



Figure 4. hi the "two-cut phase" , is zero in 
a finite neighbourhood of the origin, and shows 
two disconnected peaks. 



phase (with <I> fluctuating around zero), a uniform order phase (in which <I> typically fluctuates 
around either of the two classical minima of the potential), and the new, intermediate, non- 
uniform order phase. In particular, the critical line separating the disorder phase from the 
non-uniform order phase was identified, and the critical value for fi^ was found to scale like: 



2 0.56iV 

LI — 



In this work, it was also pointed out that the non-uniform order phase can be interpreted in 
terms of a pure potential model [871 EH], which can be obtained from equation ([3]) neglecting the 
kinetic term. The argument goes as follows: The minima of the potential term in equation ^ 
belong to orbits with representatives of the form: 



(p) 



9 



with: p G {0,1,2, ... ,iV}. 



Unless the kinetic term lifts the degeneracy associated to the p £ {1,2, . . . , N — 1} orbits, they 
are dominating - by virtue of their larger volume in the phase space - and imply non-maximal 
contributions to the expectation value of the matrix trace modulu^. 

The results display a well-behaved collapse of data when is properly scaled, and are 
(qualitatively and approximately also quantitatively) consistent with the predictions in [64] : 
the non-perturbative data seem to lie in between the theoretical expectation of [64] (which, in 
D = 2, is only an approximate one), and the pure potential model limit, in which the kinetic 
term is completely neglected. 

The results of this study were later extended in [53] , where large statistics were accumulated, 
in order to pin down the other critical line which appears in the phase diagram, namely: the one 
separating the non-uniform order from the uniform order phase. There, the phase diagram of 
the model was obtained, in terms of well-suited scaling combinations of the /x^ and g parameters 
with powers of N. In particular, the two different transition lines meet at a triple point. In 
the large- limit, the critical line between the disorder and the non-uniform order regimes is 
expressed by a relation of the form: g oc (^^^ij , whereas the transition curve separating the 



It is worth noting, however, that the dominating configurations in the model are actually those characterised 
by a smooth eigenvalue distribution. 
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Figure 5. Numerical results for the eigenvalue 
distribution in the "one-cut" regime, from data 
obtained in preliminary, low-statistics runs with 
N = 9, R^l, g = 972 and = q. 



Figure 6. As jj,^ is reduced and the transi- 
tion point is approached, the eigenvalue density 
around (f> = decreases. The data shown corre- 
spond to /i^ = —74.25, at the same N, R and g 
values as in Fig. O 



uniform order from the non-uniform order regimes appears to be compatible with a straight line 
(remarkably, the collapse of data obtained in terms of the simple scaling laws of g and /x^ works 
very well, even down to matrix sizes which are as small as N = 2). 

In the same paper, the behaviour of the specific heat per degree of freedom was worked out, 
and compared with the theoretical prediction from the large- limit of the pure potential model. 
The data for increasing matrix sizes appear to approach the expected behaviour, which would 
imply a discontinuity of the first kind in the derivative of the specific heat p.d.o.f. with respect 




A recent work [59j has refined the numerical analysis of this model, by means of a new 
method, which reduces the correlation problems affecting simpler simulation algorithms. In the 
regime, this allowed to observe an indirect signature of the non-commutative anomaly, 
as a distortion of the dispersion relation for non-vanishing momenta. In the opposite regime, 
the high-precision numerical results allowed to highlight a more detailed interpretation of the 
matrix configurations associated with I > modes. In particular, the latter are not necessarily 
incompatible with a commutative setting: non-uniform configurations would also contribute 
the path integral of QFT on the finite-radius commutative sphere, and mediate the tunneling 
among the classical potential minima. Therefore, their signature in results of simulations in the 
large- limit at fixed R is plausible, and not directly related to the UV/IR mixing. 

The situation, of course, is different, once one investigates the double scaling limit; in fact, as 
we discussed above, the meaning of the latter consists in "removing" the regularisation cut-off 
N, while keeping the non-commutativity parameter fixed. In this case, the UV/IR mixing is 
expected (as the theory tends towards an approximation of the non-commutative plane), and 
should be responsible for effects which do not have a commutative counterpart. This description 
in the double scaling limit is also compatible with the picture of the pure potential model |87^ l88j. 
as it was discussed in (52l [53] . 

The best probe to detect the non-commutative effects in the regime where the transition from 
the disorder to the non-uniform order phase takes place is probably the density of the matrix 
eigenvalues. It is possible that different observables may be more appropriate for other regions 
in the parameter space, but at the moment a general theoretical understanding of the whole 
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Figure 7. When reaches a critical value, p(0) Figure 8. The simulation results below the tran- 

vanishes; the data shown here were obtained at sition point (the plot refers to data at ~ —135) 

fi^ = —114.75. display a pattern compatible with the "two-cut 

phase" . . . 

phase diagram of the model is not consoHdated yet. Here, we shall only consider the transition 
from the disorder to the non-uniform order phase: for this region, a clear theoretical prediction 
exists [6l], which is formulated precisely in terms of the eigenvalue distribution. Technically, 
the measurement of the p{(j)) distribution appears to be simpler than other observables which 
could allow to locate the phase transition, and an algorithm like the one used in [59], which is 
highly efficient in the exploration of the phase space of the model, allows to obtain significant 
results within a limited amount of CPU time. 

In fact, the study of the eigenvalue distribution in the double-scaling limit is currently in 
progress, and the complete results will probably be published in an (upcoming) second version of 
that paper. The preliminary results obtained so far look very encouraging, showing a behaviour 
which confirms the predictions in [64j . 

As an example, in Figs. [5] to [10] we display the plots obtained from (low-statistics) inde- 
pendent data produced by test runs. These simulations correspond to = 9, holding the 
"non-commutativity parameter" = | fixed and for constant g = 972, with /j? taking four 
different, decreasing values. Note that the eigenvalues are not rescaled to the [—1,1] range, 
therefore the normalisation for the (p axis is different with respect to [M] and to Figs. [U to [H 
Also, the errorbars and complete details of the statistical analysis are not displayed for these 
preliminary data - they will be presented at the moment of the publication of the final results. 

The eigenvalue distribution plotted in Fig. [5] is obtained for /i^ = 0: its features appear to be 
compatible with the "one-cut" profile expected above the transition point (see Fig. [2]) . 

For smaller values of fi^ (Fig. the eigenvalue density around (p = decreases, and ap- 
pears to eventually vanish at a particular critical value: Fig. [7] shows the distribution profile 
approximately at the transition poinlH. 

Below the critical value, a "two-cut" pattern emerges: the numerical results shown in Fig. [8] 
seem to confirm that p{(f)) becomes peaked around two (symmetric) values, and, more important, 
that a gap of zero probability opens up for eigenvalues in a neighbourhood of zero. This is more 

^It should be remarked that, in order to identify the transition properly, one should perform a study of 
the scaling properties, when the parameters are varied; this will be addressed elsewhere. Here, we shall not 
make any quantitative claims about the precise location of the critical lines. Still, the present discussion is (at 
least qualitatively) consistent with the expected theoretical pattern, and it provides the preliminary information 
necessary to address the complete study. 
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Figure 9. ... and a zero-density gap is observed 
in a finite interval around = 0. 




Figure 10. For /.(^ values deeper in the "two-cut" 
regime (/.t^ = —351), the distribution gets more 
localised around two narrow peaks, and a zoom to 
the origin of this figure gives stronger evidence of 
vanishing eigenvalue density in a (larger) neigh- 
bourhood of = 0: a large range exists, in which 
no eigenvalues are observed. 



clearly visible in Fig. [9l which is a zoom around the origin of the previous plot. Note, however, 
that in this case we are on the borderline of the statistical precision allowed by our sample, and 
the scale of the vertical axis is already comparable with the minimal resolution of the data. 

In order to obtain more reliable information and see more precisely whether the eigenvalue 
distribution is indeed vanishing in this regime, it is instructive to extend the analysis to further, 
lower /i^ values. The result is displayed in Fig. [TOl which confirms the pattern observed in Fig. [HI 
showing a distribution which is more strongly localised around two sharp peaks; a zoom to the 
origin shows - at least within the precision of the data - an exactly vanishing eigenvalue density 
in a larger neighbourhood of </) = 0. 

Not surprisingly, the general qualitative features emerging from this picture are also shared 
by the three-dimensional model discussed in [5l], in which two non-commutative directions are 
combined with a (compactified) commutative "time" , and the regularised model is described by 
the action in equation ([7|). Indeed, also in that case, a phase was observed, which is unknown 
in the model in ordinary commutative space, and the different limits for R and N reproduce 
different physical theories. 

4 Conclusions 

To summarise the discussion presented above, we can conclude that the fuzzy space formulation 
encodes a very rich structure, describing different physical models. 

On one hand, since the existence of "anomalous" effects, which have no commutative coun- 
terpart, shows that the naive formulation of the model would not reproduce the expected com- 
mutative theory, it would be interesting to address a detailed study of some "improved" version 
of the matrix model, which may yield the correct commutative limit. 

On the other hand, the use of fuzzy spaces as a regularisation scheme for intrinsically non- 
commutative models appears to be one of the most powerful tools to investigate their non- 
perturbative features, both analytically and numerically. 
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After reviewing the main problems in fuzzy models, we presented a preliminary sample of 
new numerical results, which are in agreement with the pattern predicted by the theory for 
the phase transition between two different regimes in the scalar model on the fuzzy sphere. 
In particular, we showed evidence supporting the picture of a phase transition which can be 
investigated looking at the topological properties of the matrix eigenvalue distribution. 

Although in this paper we focused our attention onto some specific issues, the problems 
related to fuzzy spaces are among the main research interests of many different groups worldwide; 
therefore we apologise with the authors whose results have not been mentioned or properly 
discussed here. A lot of work has been done, and a lot more is currently in progress in this 
exciting research field - a field which, in view of the present status of the research, looks very 
promising. 
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